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Os ; ABSTRACT 

We describe a modification of the usual definition of astronomical magnitudes, 
replacing the usual logarithm with an inverse hyperbolic sine function; we call 
these modified magnitudes 'asinh magnitudes'. For objects detected at signal- 
to-noise ratios of greater than about five, our modified definition is essentially 
identical to the traditional one; for fainter objects (including those with a formally 
negative flux) our definition is well behaved, tending to a definite value with finite 
errors as the flux goes to zero. 
q ■ This new definition is especially useful when considering the colors of faint 

objects, as the difference of two 'asinh' magnitudes measures the usual flux ra- 
tio for bright objects, while avoiding the problems caused by dividing two very 
Oh 1 uncertain values for faint objects. 

The Sloan Digital Sky Survey (SDSS) data products will use this scheme to 
+- » ■ express all magnitudes in their catalogs. 

> 

• i— i . 

^ Subject headings: catalogs, methods: statistical, techniques: photometric 

1. Introduction 



The advantages of using a logarithmic scale to measure astronomical fluxes are obvi- 
ous: the magnitude scale is able to span a huge dynamic range, and when relative colors 
are needed they can be computed by simply differencing magnitudes measured in different 
bandpasses. These advantages are quite clear for bright objects where noise is not an issue. 

^or the SDSS collaboration. 
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On the other hand, as fluxes become comparable to the sky and instrumental noise, the cor- 
responding magnitudes are subject to large and asymmetric errors due to the singularity in 
the magnitude scale at zero flux; indeed, if a noisy measurement of an object's flux happens 
to be negative, its magnitude is a complex number! Astronomers have generally handled 
these cases by specifying detection flags and somehow encoding the negative flux. These 
problems become even more pronounced as we work in multicolor space. An object can be 
well detected and measured in several of the bands, and still fail to provide a measurable flux 
in others. In flux space the object would be represented as a multivariate Gaussian prob- 
ability distribution centered on its measured fluxes, not necessary all positive. Important 
information is lost by simply replacing a flux by an arbitrary 2- or 3- a upper limit. In magni- 
tude space, the error distribution of such an object has an infinite extent in some directions, 
making meaningful multicolor searches in a database impossible; all the non-detections in 
one or more bands must be isolated, and treated separately. One solution to this problem is 
to use linear units (e.g. Janskys), although this makes studies based on flux ratios (colors) 
inconvenient. This paper proposes an alternative solution, a modification of the definition 
of magnitudes, which preserves their advantages while avoiding their disadvantages. 



2. The Inverse Hyperbolic Sine 



We propose replacing the logarithm in the traditional definition of a magnitude with an 
inverse hyperbolic sine function. This function becomes a logarithm for large values of its 
argument, but is linear close to the origin. 



sinh 1 (x) = In x + V x 2 + 1 



sgn(x) In \2x\, if |x| 1 
x, if \x\ 1 



(1) 



The usual apparent magnitude m can be written in terms of the dimensionless normalized 
flux x = f/fo as 

m = —2.5 log 10 x = —(2.5 log 10 e)\nx = —a liax. (2) 

where / is the flux of an object with magnitude 0.0, and a = 2.5 log 10 e = 1.08574 is Pogson's 
ratio (Pogson 1856). 



Let us define the new magnitude \i as 

fi{x) = - 



sinh" 1 ^ ) 1- lub 



(3) 



Here a and b are constants, and b is an arbitrary 'softening' which determines the flux level 
at which linear behaviour sets in. After some discussion, we have adopted the name 'asinh 
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magnitudes' for /i. Consider the asymptotic behaviour of /i, for both high and low x: 
lim u(x) = — alnx = m lim u(x) = —a 

x^oo r y ' x^O 

Thus for x — > oo, /i approaches m, for any choice of b. On the other hand when \x\ <^b, 
fi is linear in x; for x <C —b, we gradually return to logarithmic behaviour, although this 
regime is never of astronomical interest for any reasonable choice of b. Intuition suggests 
that b should be chosen to be comparable to the (normalized) flux of an object with a signal- 
to-noise ratio of about one; the following section discusses the choice of b, and the related 
question of /i's error distribution. 



x 
2b 



+ lnb 



(4) 



3. The Errors in /i and the choice of the Softening Parameter b 



Although the softening parameter b can be any positive number, we show below that 
one particular selection has a couple of attractive features. In making our choice, we use 
the following guiding principles: 1) Since asinh magnitudes are being introduced to avoid 
the problems that classical magnitudes manifest in low signal-to-noise data, the differences 
between asinh and classical magnitudes should be minimized for high signal-to-noise data. 
2) We should minimize /i's variance at low flux levels. The latter is not strictly required as a 
too-small value of b merely stretches out /i's scale along with its errors, but it is convenient 
if /i's variance at zero flux is comparable to its variance at a signal-to-noise ratio of a few. 

In reality, our measurement of the normalized flux x will be noisy, with variance a 2 . We 
wish to choose b to minimize /i's variance, while keeping the difference m — ji small. Let us 
therefore compute the variances of m and fi (keeping only the linear terms in their Taylor 
series), and also their difference; arrows indicate the asymptotic behavior as x — > 0: 



Var(m) = 
Var(/t) = 

/i = a In 



a 2 a 2 



a 2 a 2 



m 



Ab 2 + 



4b 2 jx 2 




(5) 



What are the disadvantages of taking either too low or too high a value for bl Choosing 
a low value of b causes the difference m — /i to become smaller, i.e. the two magnitudes track 
each other better, but unfortunately /i's variance at x = varies as 1/b 2 . Choosing too high 
a value has the opposite effect: the difference explodes at low values of x, simply due to the 
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singularity in the logarithm in the definition of m. At the same time, /i's variance remains 
small. 

In order to balance these two competing effects, we shall determine b by minimizing a 
penalty function containing terms due to both, added in quadrature. 

The difference between the two magnitudes, normalised by ra's standard deviation, is 
given by 



m(x) — n(x) _ b fx 



^/Var(m) 



a 



5(x) = — -, = —F ( — ) , where F(y) = yln 



l + Jl + 4/n 2 



(6) 



The function F(y) has a maximum value of approximately 0.5 at y — 0.7624, so the largest 
possible deviation between the two magnitude scales is 



(m - /i) max « — ^/Var(m) (7) 

The other 'cost' associated with the choice of b is the size of the error box for fi at x = 0, 
which is 



Var^)!^ = ^ 



The total penalty can be obtained by adding these two costs in quadrature: 



<Wx + Var(^)| x=0 = — ^ + 



o 2 a 2 cr 2 a 



4a 2 46 2 4 



6 2 acr 2 



+ 



aa 2 6 2 



(9) 



which has the obvious minimum at b 2 = aa 2 . Thus the optimal setting is the value b = 
y/aa = 1.042a. As expected, b is approximately equal to the noise in the flux. This choice 



of b leads to m — fx having a maximum value of 0.52y Var(m), implying that the difference 
between the two magnitudes is always smaller than the uncertainty in m (see figure [I]). If 
the error in its measured flux is la, the error in \x is ±0.52. If the flux errors are Gaussian, 
so, to leading order, are the errors in \i as the transformation from counts to \x is linear for 
\x\ £ b. 

Figure |2] shows m and \i as a function of signal-to-noise ratio, for this choice of b, along 
with their 1-a errors. 



4. Application to Real Data 



We have been working in terms of x = f / fo, but it is usually more convenient to use 
the measured fluxes directly; In terms of the non-normalised flux /, the expressions for m, 
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/i, and Var(/i) become 

m= m - 2.51og 10 /, (10) 

// = (mo - 2-5 log 10 b') - a sinh" 1 Cf/26') (11) 

and 

VarOi) = « ^ (12) 

where m = 2.5 log 10 (/ ), = /(A an d °"' = /o " are measured in real flux units (e.g. counts). 

An object with no measured flux in a given band has a \i value of m — 2.51og 10 6' 
(equation 111) , in other words the classical magnitude of an object with a flux of b'. We note 



in passing that this value /i(0) is a convenient measure of the depth of a survey, containing 
information about both the noise properties of the sky and the image quality. 

In the discussion above we considered the idealized case of all objects having the same 
error, dominated by sky noise. This case covers most objects found in a given deep survey, 
as most are detected at the flux limit and are typically at most marginally resolved. For 
bright objects, of course, the difference between m and \x is entirely negligible. 

The optimal measure of the flux of a faint stellar object is given by convolving its image 
with the PSF. If the noise is dominated by the sky and detector, the variance of the measured 
flux is independent of the object's brightness, and is given by the background variance per 
unit area multiplied by the effective area of the PSF (47ra 2 if the PSF is Gaussian with 
FWHM 2\/2 In 2 a). If we decide upon a typical seeing quality and sky brightness for a given 
band, this defines cr's nominal value, <7o, which sets b once and for all. Each band has its 
own value of b. 

As observing conditions change so do measurement errors, with the result that the error 
in fi for very faint objects is no longer exactly the 0.52 magnitudes that it would be under 
canonical conditions. Whenever a precise error is needed for a given object's /i, it may be 



found by converting /j back to flux, or by applying equation |T^; for faint objects this reduces 
to multiplying the quoted error by g/ctq. Failure to apply such a correction would mean that 
the quoted errors on /j were wrong. 

It would be possible to choose b separately for different parts of the sky, but this would 
make the conversion of \x back to flux impractically complicated, and a significant source 
of mistakes for users of the data. The behaviour of /i as 6 changes is reasonably benign; 
the error at zero flux varies only linearily with b, as does the flux where m and /x begin to 
diverge. 

Care is also required whenever the measured flux has different noise properties, for 
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example if the flux is measured within a circular aperture or a given isophote rather than 
using a PSF. In this case, the appropriate value of a may be much larger from the one used 
to set b, with the consequence that the error in /i at zero flux considerably underestimates 
the true uncertainty (the other case, where the effective aperture is smaller than the PSF, 
is unlikely to occur in practice). It would, of course, be possible if confusing to choose a 
different set of b values for each type of (fixed size) aperture, although it seems unlikely that 
this would really be a good idea. As the discussion in section |] showed the consequences of 
even a grossly incorrect value of b are not catastrophic; the asinh magnitudes still reduce to 
our familiar magnitudes for bright objects, and are still well defined for negative values. 

One place where special care will be needed is in measures of surface brightness, where 
m and /x can depart quite strongly from one another even at levels where the flux is well 
determined. It may prove desirable to use a different value of b for such measurements; they 
are after all never directly compared with total magnitudes. 

Fan et al. (1999a) have used the asinh magnitude system to search for high-z quasars 
in preliminary SDSS data; examples of color-color plots employing asinh magnitudes may 
be found in Fan et al. (1999b). 

5. Asinh Magnitudes and Colors 

The ratio of two low signal-to-noise ratio measurements (for example, an object's color) 
is statistically badly behaved (indeed, for Gaussian distributions if the denominator has zero 
mean, the ratio follows a Cauchy distribution and accordingly has no mean, let alone a 
variance!). What is the behavior of our asinh magnitudes when used to measure colors? 

For objects detected at high signal-to-noise ratio, the difference in /i measured in two 
bands is simply a measure of the relative flux in the two bands. For faint objects this is no 
longer true, although the difference is well behaved. A non-detection in two bands has a well 
defined 'color' (/xi(0) — 7x2(0)) . As discussed above, the error in this color is approximately 
0.75<t/o"o magnitudes, assuming independent errors in the two bands. Equivalently, such a 
non-detection can be represented by an ellipsoid in multi-color space, centered at the point 
corresponding to zero flux in all bands, with principal axes 0.52ct/cto (in general a/ao will 
be different in each band). 

As an illustration of the instability of the traditional definition of color for faint objects, 
consider two objects that have almost identical colors but which are near the detection 
limit of a survey. Their asinh colors will be very similar, but (due to the singularity of the 
logarithm as the flux goes to zero) their classical magnitudes may differ by an arbitrarily 
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large amount. 

Figure |3| shows the results of a simple Monte-Carlo simulation. We took a set of 'objects' 
with 2.51 times as much flux in one band in the other (one magnitude), added Gaussian 
noise of fixed variance to each measurement, and tabulated the color measured using both 
classical and asinh magnitudes. The left hand panel shows the flux ratio, Am = mi-m2, and 
A/i = fii — fi2 as a function of signal-to-noise ratio; the right hand panels show histograms of 
their distribution in the range 1 <= S/N <= 3. At the right side of the plot, where the noise 
is less important, both Am = mi — m 2 and A/i = — /j 2 tend to -1, the correct value. As the 
noise becomes more important, the errors on the Am plot grow (and an increasing fraction 
of points in the left panel is simply omitted as their fluxes are zero or negative). The A/i 
plot shows the 'color' tending to its value at zero flux, in this case 0.0, as the signal-to-noise 
ratio drops. 



6. Summary 

We have shown that an innovative use of inverse hyperbolic sines for a new magni- 
tude scale can overcome most deficiencies of traditional magnitudes, while preserving their 
desirable features. The defining equations are 

fj, = (mo- 2.5 log 10 V) - a shffT 1 (f /2b') = /i(0) - a sinh -1 (f /2b') 

and 

Var(/i) = 4 fe' 2 +/ 2 ~ 

where a = 2.51og 10 e, / is the measured flux, a' the error in / due to the sky and detector, 
and b' a softening parameter (Equations |TI] and |T2"D . 



The principal advantages of these 'asinh' magnitudes are their equivalence to classical 
magnitudes when errors are negligible, their ability to represent formally negative fluxes, 
and their well behaved error distribution as the measured flux goes to zero. For high signal- 
to-noise ratios the difference of two asinh magnitudes is a measure of the flux ratio, while for 
noisy detections it becomes the statistically preferable flux difference; this allows meaningful 
color cuts even when an object is barely detected in some bands. Additionally, /i(0) provides 
a convenient way of summarizing the photometric depth of a survey for point-like objects, 
containing information about both the noise level of the system and the image quality. 

Asinh magnitudes will be used in the SDSS catalogs. 
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Fig. 1. — The behavior of m — /i as a function of signal-to-noise ratio, xja. The solid line 
is the value of m — /i and the region between the dotted lines corresponds to the ±la error 
region for m. The dashed lines are drawn at ±0.01. 



- 10 - 




Signol to Noise Ratio, x/o 



Fig. 2. — The behavior of m and /i, and their respective errors, as a function of signal-to- 
noise ratio, xja. The solid line is the value of \i and the region between the dotted lines 
its ±lo" error region; the points with errorbars are the classical magnitudes, m. We have 
arbitrarily chosen a zeropoint of \i = 25.0 for an object with no flux. One other feature of 
our modified magnitudes is apparent from this figure, namely that the error band on \x is 
nearly symmetrical, while the errors in m are strongly skewed at faint magnitudes. For S/N 
ratios of less than about two, m — fj, exceeds the value 0.52-y/Var(m) quoted in the main body 
of the paper; this is due to the breakdown of the linear approximation used to calculate m's 
variance. 
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Fig. 3. — The results of a Monte-Carlo simulation of the colors of a set of objects. The 
bottom panel shows the ratio of the brighter to the fainter measurement; it is clear that 
many of the low signal-to-noise points have negative flux ratios. The center panel shows the 
same points, but now in terms of the magnitude difference of the two detections. We have 
simply omitted points for which one or both fluxes were negative, although they do appear 
in the histogram. Finally, the top panel shows the difference in the asinh magnitudes. The 
'color' at zero flux was taken to be 0.0 for this simulation. 



